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ABSTRACT 



Most extrasolar planets currently known were discovered by means of an indirect method that measures the stellar wobble caused by 
the planet. We previously studied a triple system composed of a star and a nearby binary on circular coplanar orbits. We showed that 
although the effect of the binary on the star can be differentiated from the stellar wobble caused by a planet, because of observational 
limitations the two effects may often remain indistinguishable. Here, we develop a model that applies to eccentric and inclined orbits. 
We show that the binary's effect is more likely to be mistaken by planet(s) in the case of coplanar motion observed equator-on. 
Moreover, when the orbits are eccentric, the magnitude of the binary's effect may be larger than in the circular case. Additionally, an 
eccentric binary can mimic two planets with orbital periods in the ratio 2/1. However, when the star's orbit around the binary's center 
of mass has a high eccentricity and a reasonably well-constrained period, it should be easier to distinguish the binary's effect from a 
planet. 
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1. Introduction 



In iMorais & Correial (I2008I) . we studied a triple system con- 
sisting of a star perturbed by a binary system. Our aim was 
to derive an expression for the star's radial velocity in order 
to determine whether the binary's effect could be mistaken for 
that of a planet companion to the star. This question is rel- 
evant if one or even both binary components are unresolved 
which can happen for instance when these are faint but very 
common M stars. About 7 6% of nearby main sequence stars 
are M-type dLeDrewll200l . More over, in the solar neighbor- 
hood, over 50% of G and K-stars ([Duauennoy & Mavot 1 199 it 
lHalbwachs et al.ll2003l:lEggenberger et al.ll2004l) and abo ut 30% 
of M-stars (iFischer & Marcvlll992tlDeifosse et al.ll2004l) belong 
to binary or even multiple systems. 

We previously showed that the binary system's effect on the 
star is a sum of two periodic signals with very close frequen- 
cies. These signals, if detected, could be associated with two 
planets on circular orbits around the star. However, these planets 
would have very close orbits and we should be able to discard 
them as unstable systems. Moreover, we derived expressions 
for the frequencies and amplitudes of these signals, thus were 
able to identify th e hidden binary's parameters. As explained in 
Mora is & Corre ia (2008: ), our results do not agree with previous 
work bv lSchneider & Cabreral d2006l) who studied the case of an 
equal-mass binary and concluded that its effect is a single peri- 
odic signal that mimics a planet on an eccentric orbit. 

Nevertheless, in M orais & Correial (l2008h we saw that there 
are realistic situations where we can identify one of the periodic 
signals but not the other. In this case, we may mistake the bi- 
nary's effect for that of a planet on a circular orbit around the 



star. However, we can still apply our theory to compute the hid- 
den binary's parameters and then try to detect these objects. 

Our work was based on a full three-body model but as- 
sumed coplanar motion and initial circular orbits for the star and 
binary. However, general triple-star systems are likely to have 
non-coplanar motion and eccentric orbits. Therefore, in this ar- 
ticle, we extend our study to the case of three-dimensional non- 
circular motion of the star and binary system. In Sect. 2, we 
present the model and in Sect. 3 we discuss the circumstances 
that lead to the binary being mistaken by one or more planets. In 
Sect. 4, we compare the theoretical predictions with results ob- 
tained from simulations of hypothetical triple systems, we dis- 
cuss the triple system HD 188753, and we (re)analyze the exo- 
planets discovered within binary systems. Finally, in Sect. 5 we 
present our conclusions. 



2. Modeling a star perturbed by a binary system 

2.1. Equations of motion and radial velocity 

We consider the framework of the general three-body problem: 
we assume that a star with mass M* has a nearby binary system 
with masses M\ and M2 (Fig.[TJ. 

Following IMorais & Correial d2008h . we use Jacobi coordi- 
nates, which are the inter-binary distance, r^, and the distance, 
r, from the star M* to the binary's center of mass (see Fig. [TJ. 
We also assume that p = |r&|/|r| <sc 1, 

In the auxiliary frame (Fig. 1), r = (x',y',z'), where 
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x — r cos 6 , 
y' - r sin 9 , 
z' = 0, 



(1) 
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where 




Fig. 1. A star with mass M* is perturbed by a binary system with 
masses Mi and M%, The Jacobi coordinates r b and r are, respec- 
tively, the inter-binary distance and the distance of M* to the 
binary's center of mass. The unperturbed Keplerian described 
by r is in the plane (x',y'). The unperturbed Keplerian described 
by rb is defined with respect to the (x',y') plane by the angles 
Q (relative node longitude) and i (relative inclination). The ob- 
server's frame (x, y, z) with z being the line of sight is obtained 
by rotating the auxiliary frame (x',y',z') around the x axis by an 
angle / (inclination of r's orbit with respect to (x,y)). Moreover, 
S = l{r, r b ), 6 = l{x, r), and h = L(x, r b ). 

and r b = (x',y',z' h ), where 



x' b = r b (cos £1 cos(0 b - Q) - sin Q sin^ - O) cos z) . 
y'b = r b (sin O cos(0 b - Q) + cos Q sin^ - O) cos i) . 
z' h = r h sm(0 h - Q) sin i , 



(2) 



where 6 and 6 b are the orbits' true longitudes, Q is the relative 
node longitude, and i is the relative inclination. 

The transformation from the auxiliary frame {x',y',z') to the 
observer's frame (x,y,z) is 



x — x, 

y = cos/y - sin/z' , 
z = sin/y + cos/z' . 

Moreover, we define the angle S 



by 

cos S 



(3) 

l(r, r b ) which is given 
(4) 



where r is the versor of r and f b is the versor of . 

Jacobi coordinates are useful to this problem because the un- 
perturb ed motion is a composi tion of two Keplerian ellipses. As 
seen in Morais ^fe Correia! (120081) . including terms up to first or- 
der in p, we have 



Mi + M 2 , 
-G 7. (ri, + ei,0 cos Sr ■ 



h)) 



where, unless M* » M x + M 2 , 
M^ 3 



Mi + Mi 



« 1 



(5) 



(6) 



and 



(Mi+M 2 + M i ,) 
r = -G = x 



[r + e((-3/2 + 15/2 cos 2 5)r - 3 cos 5r & )] , 



(7) 



e = 



MiM 2 2 -, 

-TP 2 < — « 1 . 



(Mi + M 2 ) 2r 4 
Thus we can write r b = r b o + e b r b i , where 

fbo = , r b0 , 



and r ■■ 



ro 



'60 

ro + eri 



where 



XMi +M 2 + M i ,) i 



-ro 



(8) 



(9) 



(10) 



Therefore the 0th order solution, r^o, is a Keplerian ellipse with 
constant semi-major axis, a b , and frequency 



n h = 



i 



G(Mi + M 2 ) 



(11) 



and the 0th order solution, ro, is a Keplerian ellipse with constant 
semi-major axis, a, and frequency 



G(Mi + M 2 + M*) 



(12) 



Furthermore, from Eqs. dTTb and ( TT2b we see that in general n <K 
n b since a b <K a (this is true unless M* » Mi + M 2 ). 

Now, the distance of the star to the triple system's center of 
mass is 



Mi + M 2 



■ r , 



Mi + M 2 + Mi, 

hence from Eq. (0 we have 

Mi+M 2 
r+ = -G , x 



(13) 



[r + e((-3/2+ 15/2 cos 2 S )f - 3 cos S r b )] 



(14) 



Since r 2 = r 2 (l + 0(e)), r = r (l + 0(e)) and f b = r M) (l + 
0(e b )), we obtain an approximation to the previous equation that 
is accurate up to 0(e), i.e. 



Mi + M 2 „ 
-G r 

r L 

Mi +M 2 
-G z e 



(15) 



3 15 

"2 + T cos ' r ° 



3 cos S r^o 



whose solution, = r*o + er+i, is thus a combination of a 
Keplerian ellipse (with frequency ri), r*o, and a perturbation 
term 



>Vi = -G 



Mi + M 2 



3 

— + 

2 



15 



cos S ro - 3 cos S f b D 



(16) 



The radial velocity is the projection of the star's barycentric 
velocity along the line of sight (defined as the z-axis in Fig. [TJ, 
i.e., V r = z*. Therefore, we only need to compute the z compo- 
nent of the first order correction, er*i (Eq. (fT6ll). i.e. 

Mi + M 2 
z* = -G z ex 



- + 

4 



'o 
15 



cos(25) ro • z - 3 cos Sr b0 -z 



(17) 
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Additionally, from Eqs. ([TJ, (0, ([3), and © we have 

r -z = sin(0)sin7, (18) 
f b0 -z = sin(0(,) sin 7 + sm(9 b - £2) x (19) 

(sin z cos 7 - (1 - cos z) cos £2 sin 7) 
cos S = f -h~h - ho 

= cos(6 b - 0) - (1 - cos /) sin(0 - £2) sin(0 fi - Q) . (20) 

2.2. Star and binary on circular orbits 

In this case , ro - a, r b Q — a b , — A — nt + Aq, 6 b - A b - 
n b t + Aw, where A and A b are the orbits' mean longitudes. Thus 
Eq. ( fTTT i becomes 

z* = - 5n\ sin 7 x [4( 1 - 3 cos/ 2 ) sin(« f + Aq) 

+2(1 - cos i 1 ) sin(n f + A - 2 £2) 

-10(1 - cos i 2 ) sin(3 n f + 3 Aq - 2 £2) 

+5(1 + cos z) 2 sin((2 nj, - 3 «) ? + 2 Am ~ 3 Ao) 

-5(1 - cos 2 ) sin((2 n b + 3 n) t + 2 A w + 3 /lo - 4 Q) 

-(1 + cos z') 2 sin((2 n b - n) t + 2 A b Q - Aq) 

+6(1 - cos z' 2 ) sin((2 n b - n) t + 2 A b Q — Aq — 2 £2) 

+(1 - cos i) 2 sin((2 n b + n) t + 2 Ajo + Ao - 4 £2) 

-6(1 - cos z' 2 ) sin(2 n b + n) t + 2 Am + Aq - 2 £2)] 

+6 <5 n 2 b cos 7 x [2 sin z cos i sin(« f + Aq - £2) 

+ sin z'(l + cos z) sin((2 n b — n) t + 2 A^o - Ao — £2) 

+ sin z'(l - cos z) sin((2«i + n)t + 2 At® + Ao - 3 Q)] ,(21) 

where 

6= 8(M 1+ M 2 )47J Qfe (22) 

as defined in M orais &Correiald2008l) . 

The right-hand side of Eq. ( |2TI ) is a sum of periodic terms, 
hence V r = z* is obtained by integrating them with respect to 
time. Therefore, z* is a linear combination of six periodic terms 
with frequencies n,3>n,2n b -n, 2n b + n,2n b -3 n, and 2 n b + 3 n. 

In the coplanar case (z = 0), we recover the solution obtained 
in lMorais & Correial d2008l) 

n 2 

z* = 66 sin 7 — cos(nf + Ao) 
n 

n\ 

-15(5 sin 7- — cos((2n^ - 3 n) t + 2 Aq - 3 Ami) 

2 n b — 3 h 

«2 

+36 sin 7- — 2 — C os((2 n b - n) t + 2 A - A M ) (23) 
2«i — n 

i.e., is a composition of three periodic terms with frequencies 
n, 2 n b — n, and 2 n b — 3 n. 

We note that when z = (coplanar case) the angle 7 defines 
the projection of the star's orbit along the line of sight and thus 
z (and z) scale with sin 7. However, when i + the angle that 
defines the projection of the star's orbit along the line of sight 
depends on 7, z, and £2. When z + 0, the amplitudes associated 
with the frequencies n and 2n b + n indeed do not scale with sin 7 
(cf. last three terms in Eq. (OTt ). In particular, when 7 = (pole- 
on configuration) and if i + we can still detect these terms 
(frequencies n and 2n b + n), which are due to the binary's effect 
on the star. 
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Fig. 2. Circular case with 7 = 90° (equator-on configuration). 
Normalized amplitudes A of periodic terms as function of the 
relative inclination z. These are obtained by integrating Eq. (OTT i 
with respect to time. The true amplitudes are A6n b . 



In Fig. (2), we show the normalized amplitudes (A) associ- 
ated with the frequencies 2 n b + n and 2n b ± 3 n, as a function of 
z for £2 = 0, 7 = 90° (Fig. 2, top) and £2 = 90°, 7 = 90° (Fig. 2, 
bottom). In Fig. (3), we show the equivalent picture for 7 = 0. 
The true amplitudes are A6n b . 

At 7 = 90° (equator-on configuration) and small z, the term of 
frequency 2 n b — 3 n is dominant but as z increases its amplitude 
decreases, while the amplitudes of the terms with frequencies 

2 n b + 3 n and 2n b ±n increase until they are all approximately 
equal at i = 90°. At 7 = (pole-on configuration), only the 
terms with frequencies n and 2n b ±n appear (cf. Eq. (I2TI) ). hence 
we cannot detect the star's motion around the binary's center of 
mass but if i + we can still detect the binary's effect on the star. 
Moreover, the amplitudes associated with the frequencies 2n b + 

3 n are independent of £2 at any value of 7 while the amplitudes 
associated with the frequencies 2 n b ± n depend on the orbits' 
intersection £2 except when 7 = (cf. Eq. (f2TJ). The pictures for 
90° > i > 0° (star and binary on prograde orbits) and 180° > z > 
90° (binary on a retrograde orbit) are symmetrical except that the 
frequency n is replaced with -n, which corresponds to inverting 
the star's (or binary's) motion. 

2.3. Star and binary on eccentric orbits 

We assume that the Keplerian orbit, ro, has a semi-major axis 
cz, eccentricity e, and true longitude 6 = / + m (where / is the 
true anomaly and m is the longitude of periastron) and that the 
Keplerian orbit, r^o, has semi-major axis a b , eccentricity e b , and 



4 
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circular, I = 




60 80 100 120 
relative inclination (i) 



Fig. 3. Circular case with 1 = (pole-on configuration). 
Normalized amplitudes A of periodic terms as function of the 
relative inclination i. These are obtained by integrating Eq. ( TZTT i 
with respect to time. The true amplitudes are A 6 rib. 



true longitude 6b = fb + tub (where fb is the true anomaly and 
mb is the longitude of periastron). Thus Eq. (fTTl i becomes 



z* = -G- 



a- 
15 



M 2 I n 
— eo 

" r 



— ) x 

ab ) 



cos(2S ) \h ■ z - 3 cos S ho • z. 



where 



eo = 



M\M 2 
(M, + M 2 ) 2 



(24) 



(25) 



We can express the orbital solutions, ro and r^o, as Fourier 
series of the mean anomalies M = n t + Mq and Mb = rti,t + Mi,o, 
respectively, also known as elliptic expansions. In particular, we 
obtain expressions for (q/rp ) 4 and (rbo/ab) 2 using ellip tic expan- 
sions for r/a and rb/cib (e.g. lMurrav & Dermotll(ll999l) \ We ob- 
tain similar expressions for % ■ z (Eq. (fTSTl). ho • z (Eq. (fT9l), 
and cos S (Eq. ( l20l i) usi ng elhptic expansions for c os(/), sin(/), 
cos(/i), and sin(/ fc ) (e.g. lMurrav & Dermottl (1 19991) 1. 

2.3.1. Coplanar motion 

In this case, i = 0, thus from Eqs. (fl8t . ( fT9] l, and ( l20t 



ro • z = sin(#) sin / , 
h-Z = sin(Gb) sin / , 
cos S = cos(8i, - 6) , 



(26) 
(27) 
(28) 



where 6 — f + m and 6b — fb + ?Vb- 

The next step is to replace these expressions into Eq. (l24l 
using elliptic expansions for (a/ro) 4 , (rbo/ab) 2 , sin(6>), sin(0/,), 
cos(6b - 6), and cos(2 (6b - 6)). We start by using elliptic expan- 
sions up t o first order in the eccent ricities. 

As in iMorais & Correia (l2008l) . our aim is to obtain a solu- 
tion to Eq. (l24T i that is a linear combination of periodic terms. 
However, when e + we have a constant term of the type 
esin(cr) appearing in z\- We can remove this secular term by 
rewriting Eq. ([ToT l as 



^M\+M 2 



-p + r *! 



(29) 



where the vector p = (p x ,p y ,p z ) has radius (3/4) e, and is 
aligned with the periapse of the unperturbed orbit, ro, i.e. 



Px - (3/4)ecos(nr) , 

p y = (3/4)esin(tir)cos/ , 

p z = (3/4)esin(cr)sin/ , 

and 

Mi+M 2 Mi+M 2 
r +i = G p-G X 



^ + — cos(25) ) i\> - 3 COS.V i 



(30) 
(31) 
(32) 



(33) 



Now, the first term in Eq. d29| i can be interpreted as the slow 



1 II 

rotation (frequency ~ eJ n <K n «: «/,) of the unperturbed or- 
bit ro's periapse described by the vector p. Therefore, we can 
concentrate in obtaining the solution to Eq. d33l . 

The z component of the first order correction, er'*i given by 
Eq. d33l ), becomes 

z* = 3 6 n\ sin / x [-2 sin (n t + A$) 
— sin ((2 rib - n) t + 2 ^bo - M)) 
+5 sin ((2 rib - 3 n) t + 2 Abo - 3 Ao) 
-15 eb sin((«2, — 3n)t+ Abo - 3 Aq + vib) 
-e b sin ((3 n h - n) t + 3 A b o - A - m b ) 
+5 eb sin ((3 nt, — 3 n) t + 3 Abo — 3 A® - mb) 
+2 e b sin {(rib + n)t + Abo + Ao - mb) 
+3 e b sin ((n h -n)t + A b o~ Aq + m b ) 
-2 e b sin ((n h -n)t + A b o- A - m b ) 
-6 e sin (2 n t + 2 Aq - w) 
—e sin (2ribt + 2 Abo — w) 
+25 e sin ((2 «/, - 4 «) f + 2 Abo - 4 Aq + w) 
-3 e sin ((2 rib - 2 n) t + 2 Abo — 2Aq + m) 
-5 e sin ((2 rib - 2 n) t + 2 Abo -2Aq — vj)] , (34) 

where Aq - M + m and Abo = Mbo + ^b- 

The right-hand side of Eq. (t34b is a sum of periodic terms, 
hence V r = z* is obtained by integrating these with respect to 
time. Therefore, when including terms up to first order in the 
eccentricities, z* is a composition of 12 periodic terms with fre- 
quencies: n,2rib-n and 2 rib - 3 n (zero order); rib - 3 n, 3 rib - n, 
3 rib — 3 n, rib + n and rib ~ « (first order: e^,); 2n,2rib,2rib — 4-n 
and 2 - 2 n (first order: e). The amplitudes of the terms with 
frequencies rib - n and 2rii, - 2n depend on zdt, and cr, respec- 
tively. 

We can easily extend our theory by including higher-order 
eccentricity terms in the elliptic expansions[3(i.e., Fourier series 
of the Keplerian orbital solutions ro and r^o). In particular, if we 
include terms up to second order in the eccentricities, there are 
additional periodic terms with frequencies: 3 n, 2 rib + n, 4 ri/, - n, 
4rib-3n (second order: e?); 3 n,2rib+n, 2n b -5n (second order: 
e 2 ); 3rib — 2n, rib, 3nb,nb — 4n,nb + 2n,nb-2n,3rib — 4n (sec- 
ond order: eeb). Moreover, the amplitudes associated with the 
frequencies 2ni, - n and 2 rib - 3 n (zero order) have corrections 
of second order in the eccentricities. 

In Fig. (4), we show the largest normalized amplitudes (A > 
2) of the periodic terms that appear up to fourth order in the 



1 We performed these expansions using the computer algebra soft- 
ware Maple. 
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Fig. 4. Coplanar case (i — 0) with star on an eccentric orbit 
(e = 0.1). Normalized amplitudes A of periodic terms, as func- 
tion of the binary's eccentricity e/,. We show all periodic terms 
that appear up to fourth order in the eccentricities and reach at 
least A — 2. The true amplitudes are A6n\, sin /. 



Fig. 5. Coplanar case (i = 0) with star on an eccentric orbit (e = 
0.1). Normalized amplitudes A of periodic terms, as function of 
the star's eccentricity e. We show all periodic terms that appear 
up to 10th order in the eccentricities and reach at least A — 2. 
The true amplitudes are A 6 tit, sin /. 



eccentricities (except harmonics of n) as functions of e/, when 
e = 0.1. The expansion in eb converges rapidly thus it is not 
necessary to include higher order terms at least up to e\, = 0.4. 
When eb is small, the term with frequency 2nj — 3 n is dominant 
but when eb > 0.16, the term with frequency rib — 3n has the 
largest amplitude. Other important terms have frequencies 2 rib - 
An and rib - An. 

In Fig. (5), we show the largest normalized amplitudes (A > 
2) of the periodic terms that appear up to 10th order in the eccen- 
tricities (except harmonics of n) as functions of e when e/, = 0. 
We include corrections in the amplitudes up to 12th order in the 
eccentricities. These high order expansions are necessary to ob- 
tain results valid up to e — 0.4 since the expansion in e converges 
slowly. When e is small, the term with frequency 2 rib ~ 3 n has 
the largest amplitude but as e increases, the dominant term be- 
comes 2 nb - A n (e > 0.18), 2 rib — 5 n (e > 0.28), and 2 «/, - 6 n 
(<? > 0.36). As e increases, the number of terms with frequencies 
near 2 rib an d similar amplitudes increases. 

2.3.2. Non-coplanar motion 

When i + 0, we must follow the same procedure described in 
the previous section, expressin g the orbital solutions rp a nd r^o 
using elliptic expansions (e.s. iMurrav & Dermottl dl999t) ">. but 
now using the general expressions for ro • z, ?bo • z, and cos S 
(Eqs. dTSj, (2D, and ©). 

We can show that, including terms up to first order in the 
eccentricities, z* (and z*) are linear combinations of 21 periodic 
terms with frequencies n, 2 n, 3 n, A n, 2 nb, 2 nb+n, 2 n^+3 n, nb+ 
n, rib±3 n, 3 nb+n, 3 n/,+3 n, 2 n^+4 n, and 2 n/,+2 n. We note that, 
except for the additional harmonics of n, these include the set of 
frequencies already present in the coplanar eccentric case plus 
an additional set of frequencies obtained from the previous set 
by replacing n with -n, which as noted previously corresponds 
to inverting the star's (or binary's) motion. 




60 80 100 120 
relative inclination (i) 



Fig. 6. Eccentric case (eb = 0.1) with / = 90°. Normalized am- 
plitudes A of periodic terms as a function of the relative inclina- 
tion i. The true amplitudes are A 6 rib- 



We already described how the terms with frequencies 2nb + 
n, and 2«^ ± 3« (circular case) behaved with /, i, and Q. 
(Figs. (2), (3) and Sect. 2.2). To understand what happens when 
e and eb are small but non-zero, we now describe the behavior 
of the terms that appear to first order in the eccentricities. 

When eb + 0, the frequencies rib ± 3 n, 3 rib ± n, nb ± n, and 
3nb + 3n appear. In Fig. (6), we show the normalized amplitudes 
(A) as functions of i for / = 90° (equator-on configuration), eb = 
0.1, and Q. = 90° . At / = 90° and small i, the term with frequency 
rib — 3 n has maximum amplitude. As the relative inclination i 
increases, this amplitude decreases and at i = 90° the terms with 
frequencies nb + n have maximum amplitude. 

When e + 0, the frequencies 2 n^ 2 nb ± 2 n, and 2 rib ± 4 n 
appear. In Fig. (7), we show the normalized amplitudes (A) as 
functions of i for / = 90° (equator-on configuration), e = 0.1, 
and Q. = 90°. At / = 90° and small i, the term with frequency 
2 rib - 4 n has maximum amplitude. As the relative inclination i 
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Fig. 7. Eccentric case (e = 0.1) with 7 = 90°. Normalized ampli- 
tudes A of periodic terms as a function of the relative inclination 
i. The true amplitudes are A 6 rib. 



Fig. 9. Eccentric case (e = 0.1) with 7 = 0. Normalized ampli- 
tudes A of periodic terms as a function of the relative inclination 
i. The true amplitudes are A 6 rib. 
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Fig. 8. Eccentric case (fib =0.1) with 7 = 0. Normalized ampli- 
tudes A of periodic terms as a function of the relative inclination 
i. The true amplitudes are A 5 rib- 



increases, this amplitude decreases and at i = 90° it matches the 
amplitudes of the terms with frequency 2 rib ±2 n. 

In Figs. (8) and (9), we show the normalized amplitude (A) 
when I - Q, eb =0.1, and e = 0.1, respectively. In this config- 
uration (7 = 0), only the frequencies 3 rib + n, rib + n, 2 n/„ and 
2 rib ± 2 n have non-zero amplitude. As noted previously, in this 
case (pole-on configuration) we cannot detect the star's motion 
around the binary's center of mass but if i + 0, we can still detect 
the binary's effect on the star. 

Moreover, when i + only the amplitudes associated with 
the frequencies rib + 3 n, 3 rib + n, rib + n, 2 fib, and 2 rib ± 2 n 
depend on the relative node longitude, Q, since the other fre- 
quencies originate from the angles 2 6b + 36 and 3 6 whose am- 
plitudes are independent of D. (cf. third, fourth and fifth terms in 
Eq. (fZTJ). We can also show that the amplitudes associated with 
the frequencies rib ±3n and ri\, ± n depend on the binary's longi- 
tude of pericenter Wb, while the amplitudes associated with the 
frequencies 2 rib and 2 rib ± 2 n depend on the star's longitude of 
pericenter m. 

The pictures for 90° > i > 0° (star and binary on prograde or- 
bits) and 180° > i > 90° (binary on a retrograde orbit) are sym- 
metrical except that the frequency n is replaced with -n, which 
corresponds to inverting the star's (or binary's) motion. 



3. Can a binary system mimic a planet? 

As in M orais & Correial d2008l) . we wish to identify the circum- 
stances in which we can we mistake the binary's effect for that of 
a planet, when we do not know whether a binary is present (one 
or even both components may be unresolved). We have seen that 
the star's radial velocity will have a slow Keplerian term, which 
can be associated with a mass, Mi + M 2 , located at the binary's 
center of mass. Hence, from the radial velocity we can infer the 
presence of a nearby "star", Mi + M 2 , except when 7 = (pole- 
on configuration) where we can nevertheless detect the wobble 
due to the binary if i + 0. Moreover, the radial velocity of the 
star, M*, that has a planet companion, M p , in a circular orbit 
with frequency n p , semi-major axis a p and inclination I p , and 
another nearby "star", Mi + M2, is 



V r = V + 



M p sin I p 



M* + Mi + M 2 + M„ 



n p a p cos(n„f + m p ) 



(35) 



where Vq includes the Keplerian motion around the "star" Mi + 
M2, and we can neglect M p in the denominator of Eq. d35l l since 
the planet's mass is negligible with respect to the stars. 

We saw above that the radial velocity of a star perturbed by 
a binary system is a composition of several periodic signals. The 
terms that include harmonics of the frequency n can in principle^ 
be associated with the star's slow Keplerian motion around the 
binary's center of mass (i.e., Vq). Therefore, we have 



V r = Vo+2^ k p cos(n p t + m p ) , 
where 



K p * A5n b - A 



VG(Mi + Mi) M X M 2 
8 {Mi + M 2 ) : 



-1/2 



(36) 



(37) 



and the factor A is a normalized amplitude that, as explained 
previously, can be obtained for each frequency n p . 

Each of the terms with amplitude K p could be identified 
as a planet on a circular orbit around the star with frequency 
rip - prib + qn, where p > and q > are integers. However, 
since we have a limit on the observations' precision, we can 
sometimes detect only a few periodic terms with the largest am- 
plitudes. On the other hand, since rib ^> n, terms of equal p have 



2 This is not always true since the binary's contribution may prevent 
a Keplerian fit, in particular if the data set covers the long-period T = 
2n/n. 
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very close frequencies. Hence, if we can detect these latter terms, 
we associate them with planets on very close orbits that can be 
discarded as unstable configurations. However, since we have a 
limit on the observations' resolution, we are sometimes unable 
to distinguish terms with very close frequencies. 

The radial velocity method of extrasolar planet's detection 
consists of measuring tiny Doppler shifts in the star's spectral 
lines caused by the star's motion. These Doppler shifts are col- 
lected over an observation timespan, t b s , in order to obtain the 
star's radial velocity curve. The instrument's precision is the 
smallest Doppler shift that can be detected, hence the smallest 
detectable amplitude in the radial velocity's Fourier transform. 
On the other hand, the data's resolution is the smallest frequency 
difference that can be detected, i.e., 5f — l/t b s - Therefore, to 
distinguish two terms with frequencies that differ by A q ■ n, we 
must have 



tabs > 



(38) 



Aq 

where T = 2n/n, 

Now, if all detected frequencies are well separated (i.e., have 
different p) we may then be led to believe that the star has one 
or more planet companions. In that case, we apply Kepler's third 
law to derive the fake planets' semi-major axis, a p , from the sig- 
nal's frequencies, n p , i.e. 

(GM») 1/3 

a P = 7 ^2/3 ' (39) 

and from Eqs. d35l l and (l36b we obtain the associated minimum 
masses, M p sinl p , from the signal's amplitudes K p (Eq. (l37li). 
i.e. 



M p sin l p 



K„ 



(40) 

M* + Mi + M 2 a p n p 

Moreover, from Eq. ( 1371 ) we see that the binary's effect is more 
pronounced (i.e K p is larger) for a large ratio a^ja, small inter- 
binary distance ab, and a massive binary of equal stars {M\ = 
M 2 ). 

3. 1 . Circular coplanar orbits 

This was studied in M orais & Correial d2008l) . In the case of 
coplanar circular prograde orbits, we recall that the periodic 
terms due to the binary's effect have frequencies n p = 2 rib-n and 
n p = 2rib-3n, and associated amplitudes K p » (3 /2)rib5 sin / 
and K p « (l5/2)ribSsinI, respectively. As noted previously, 
these signals may at first sight be identified as two planets or- 
biting the star with very close frequencies. However, if we could 
resolve these close frequencies (i.e., if t„b s > T/2, cf. Eq. d38])), 
we should be able to discard the planet's existence since these 
systems are most certainly unstable. There are nevertheless situ- 
ations where we can mistake the effect of a binary system for a 
planet: 

- We cannot resolve the two close frequencies (i.e., if t \, s < 
7 '/2, cf. Eq. (l38l). hence we detect only one signal at ~ 2 rib. 

- Owing to limited instrument's precision, we detect the signal 
at 2 rib ~ 3 n but not the one at 2 rib — n (these have amplitudes 
in the ratio 5/1.) 

In the latter scenario, from Eq. d39l >, we obtain the associated 
fake planet's semi-major axis 



(GM*) 1/3 
2p (2«i-3n) 2 / 3 



Mi, 



4(Mj + M 2 ) 



1/3 



and from Eqs. ( l37l l and ( l40b , since the motion is coplanar (I p 
I), we obtain the fake planet's mass 



M„ 



K„ 



15 M X M 2 

M* + Mi + M 2 ~ a„(2 n b -3n) ~ 32 (M l + M 2 ) 



a h 



.(42) 



at , 



(41) 



In the former scenario, the signal's amplitude can vary by 20% 
depending on the unresolved signals' phases. Nevertheless, these 
formulae are still approximately valid. 

3.2. Non-circular, non-coplanar orbits 

The case of inclined and eccentric orbits is similar. We saw that, 
when including terms up to first order in the eccentricities, the 
star's radial velocity is a composition of 21 periodic terms. Four 
terms are harmonics of the frequency n, hence can in principle be 
discarded as associated with the slow Keplerian motion around 
the binary's center of mass. Four terms have frequencies close 
to rib, nine terms have frequencies close to 2 rib, and four terms 
have frequencies close to 3 rib. 

Comparing Figs. 2, 3, 6, 7, 8, and 9 we see that the max- 
imum radial velocity variations associated with dominant peri- 
odic terms occur when / = 90° (equator-on configuration) and 
the orbits are nearly coplanar (z * or i « 180°). In the case 
of coplanar prograde motion (i = 0), the dominant term has fre- 
quency 2rib - 3n (for small e and e/,', cf. Figs. 4 and 5) and 
frequency «/, - 3 n (for large eb', cf. Fig. 4). In particular, when 
the binary's orbit is eccentric it is possible to detect signals at 
2 rib -3 n and rib ~ 3 ft, which can mimic two planets with orbital 
periods in the ratio 2/1. However, as e increases the number of 
terms with frequencies close to 2 rib and similar amplitudes in- 
creases (Fig. 5), hence it should be easier to detect signals with 
very close frequencies. 



4. Examples 

4.1. Theory versus simulations 

To test our model, we performed some numerical simulations 
of a triple system composed of main sequence stars with dif- 
ferent spectral types G, K, and M, and masses M* = M Q , 
Mi = 0.70 M Q , and M 2 = 0.35 M e , respectively. The smaller K 
and M stars form a binary system with semi-major axis ab = 1.1 
AU, eccentricity eb, and inclination //,. The G star is in a wide 
orbit with semi-major axis a = 10 AU, eccentricity e, and incli- 
nation / around the binary's center of mass. The M star is much 
fainter than the K or G stars, hence it represents the unresolved 
component of the binary. 

We numerically integrated this triple system and computed 
the radial velocity of the G star (the brightest body in the sys- 
tem) assuming an equator-on configuration (/ = 90°). From this, 
we simulated 100 observational data points for a time span of 
4000 days (about 1 1 years). We assumed that the data were ob- 
tained with a precision limit of 0.535 m/s, which corresponds to 
A = 4.46 (cf. Eq. [37t and is about the highest precisi on that can 
presently be obtained (HARPS, iMavor et ail d2003l) ). We then 
applied the traditional techniques used to search for planets. 

Since t b s = 4000 days and the output step is A t « 40 days, 
the largest detectable frequency is f c = 2/(A?) = 1/20 days -1 
and the frequency resolution is A / = 1 /4000 days -1 . We present 
the theoretical predictions versus simulations in Table 1 . 

In a first step, we can only detect the large-amplitude radial 
velocity variations due to the star's slow motion (with period 
T = 2n/n) around the binary's center of mass. Once we subtract 
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these long-term variations we are able to detect the signal due to 
the binary's motion, which could be mistaken for one or more 
planet companions to the star. 

Examp le 1 is a circular coplana r system (i = 0) already pre- 
sented in iMorais & Correial d2008l) . where the binary mimics a 
planet of approximately 20 Earth masses. Examples 2 and 3 are 
also circular but inclined systems with i = 30° and i = 60°, 
respectively. These demonstrate that increasing the relative in- 
clination i leads to a decrease in the planet's mass M p , as we 
would expect from Eq. ( |2"TT i since the amplitude associated with 
the frequency 2 rib - 3 n is proportional to (1 + cos z') 2 . 

The next four examples are coplanar systems (z = 0) with 
e — 0.1 and et =0.1, 0.2, 0.3, and 0.4, respectively. When et, = 

0. 1 (Ex. 4), we detect the signals with frequencies 2 rib — 3 n and 
n b -3n (cf. Fig. 4 with A = 4.46). When e b > 0.2 (Exs. 5, 6 
and 7), the signal with frequency nt, - 4 n is above the detection 
limit (cf. Fig. 4 with A = 4.46) but we cannot distinguish it from 
the signal with frequency rib-3n since the resolution is not high 
enough. To resolve a frequency difference of n we would indeed 
require a resolution of at least Af = 2n/n, which corresponds 
to a time span T = 2n/n, i.e., about 22 years. At present, such a 
large timespan is not realistic since exoplanet detections started 
only 15 years ago. This lack of resolution prevents the correct 
detection of the signal with frequency rib ~ 3 n in the 5th, 6th, 
and 7th examples since this is "contaminated" by the close signal 
with frequency rib - 4 n. 

In examples 1, 2, and 3, the binary mimics a planet while in 
examples 4, 5, 6, and 7, it mimics two planets with orbital peri- 
ods in the ratio 2/1 . In these cases, the data's precision and ob- 
servation timespan do not allow us to distinguish the binary's ef- 
fect from planet(s). We now present an example where we detect 
signals with very close frequencies thus can reject the planet(s) 
hypothesis. 

Example 8 is a coplanar system (z = 0) with e b = 0.05, and 
e = 0.4. According to Fig. 5, when e = 0.4 the terms with fre- 
quencies 2 ri{, — 6 n, 2 rib — 5 n, 2 rib — 1 n, 2 n b — 8 n, and 2 rib — 9 n 
are all above the detection level (A = 4.46). However, we only 
detect signals at approximately 2 rib- 6 n and 2 rib - 8 n as the ob- 
servation timespan, t„b s ~ T/2, is not long enough to resolve the 
other frequencies. These nearby unresolved frequencies prevent 
the correct detection of the signals at 2 rib - 6 n and 2 rib — 8 n. 
Additional Fourier analysis of the residuals show the presence of 
a signal at 2 rib - 10 n but the fitting to the signals at 2 rib- 6 n and 
2 rib — 8 n deteriorates when including this frequency. The residu- 
als are 1.7 m/s, i.e., above the precision limit (~ 0.5 m/s), which 
indicates the presence of unresolved / undetected frequencies. 

Finally, we note what happens when e is large but the ob- 
servation time span is short. We performed a new analysis of 
example 8 assuming t () b s = 2000 days and t„b s - 1000 days, 

1. e., about T/4 and T/8, respectively. In the first case, we de- 
tected signals at 243 days (2 rib - 6 n) and 28 1 days (2 rib - 1 1 «), 
with amplitudes 5.0 m/s and 2.7 m/s, respectively. However, the 
residuals after the fit are still 1.3 m/s. In the second case, we de- 
tect only one signal at 257 days (~ 2 rib) with an amplitude 6.7 
m/s,, which corresponds to a 0.3 My planet at 0.79 AU. Since 
tabs T, we cannot resolve the signals at nearby frequencies 
although these affect the measured amplitude. The residuals af- 
ter the fit are ~ 0.4 m/s, hence there is no hint of unresolved / 
undetected frequencies. 

4.2. Scaling of the observables 

In the previous section, we presented examples of several config- 
urations of a triple star system where the binary is at a = 10 AU 



from the target star. However, if the binary's stars are not both 
faint, it may be difficult to obtain precise radial velocity measure- 
ments because the target star's spectrum is likely to be contam- 
inated by light received from the binary's stars. This will cause 
spectral line blending and it may be very difficult to separate 
the target star's motion from the binary's motion. Therefore, a 
realistic scenario would be a distant binary system (e.g. 50 or 
100 AU) or a close binary (e.g. 10 AU) composed of faint stars 
(e.g. M-type). 

In Table 2, we show some examples of triple systems com- 
posed of a solar mass star and a binary system on circular copla- 
nar orbits. In examples 1 and 2, the planets could b e detected by 
the cu rrent most precise instruments (e.g. HARPS, iMavor et all 
(2003)). On the other hand, the effect of the binaries in ex- 
amples 3 and 5 is slightly below the current detection limit 
but the long periods imply that we would need a long ob- 
servation timespan to constraint the orbits. Finally, the plan- 
ets mimicked by the binaries in examples 4 and 6 are not 
detectable by the current techniques but could be accessible 
to the planned exoplanet search program EXPRESS0/CODEX 
dD'Odorico & CODEX/ESPRESSO Teamll2007l) . 

Our results can be easily scaled: 

- If we increase a and ab by a factor y, then the ratio cib/a does 
not change, K p decreases by y~ 1 / 2 (Eq.[37l). a p increases by a 
factory (Eq.[39l, and M p remains unchanged (Eq.l40b. while 
the periods increase by y 3 ^ 2 . 

- If we decrease Mi and M2 by a factor y, then K p decreases by 
y 1 / 2 (Eq. [37j, T p increases by y~'^ 2 , a p increases by a factor 
y~ 1/3 (Eq.|3"9t. and M p decreases by y 1/3 (Eq.gUJl, while T 
remains unchanged. 

4.3. The triple system HD 188753 

The system HD 188753 is composed of a star (A) and a close 
binary (Ba+Bb) with an orbital period of 155 days that orbits 
the star (A) with a 25.7 year period dKonackil 120051) . A hot 
Jupiter with a 3.35 day perio d orbit around the A star was an- 
nounced by iKonackil (120051) but this was later challenged by 
lEggenberger et alj d2007l) who claimed they could not find ev- 
idence for this planet. This shows the difficulty in identifying 
planets of stars with a nearby binary system. As noted previ- 
ously, when observing the target star to obtain its spectrum we 
are likely to also detect light from the nearby stars. In the case 
of HD 188753, the contribution of Bb is modest since this star 
is faint but we cannot ignore the contributi on of Ba, which is 
blend ed within the target star A's spectrum (Eggenbergeret al. 
[2007h . 

The parameters of this triple system are M* = 1.06 M , 
Mi = 0.96 M , M 2 = 0.67 M , a = 12.3 AU a b = 0.67 AU, 
e = 0.5, e b = 0.1, and / = 34° dKonackill200lh . We can apply 
our theory to estimate the effect of the binary on the star assum- 
ing 50 data points and t b s = 500 days, which are app roximately 
the parameters reported in Eggenberge r et al.1 (l2007h . Assuming 
that the triple system is coplanar, we obtain radial velocity oscil- 
lations with an 85 day period (frequency ~ 2 rib) and amplitude 
0.5 m/s (i.e., 5.5 times the value obtained with a circular copla- 
nar model). This can be interpreted as a fake planet of 4 Earth 
masses located at 0.38 AU. As t a b s <s T = 25.7 years, we can- 
not resolve additional frequencies near 2 rib- However, the pre- 
dictions made here will change for different values of the angle 
variables (phases and relative inclination) and if we increase the 
observation time span. Nevertheless, the precision obtained for 
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Table 1. Parameters of planets mimicked by several triple system configurations when the instrument's precision is 0.535 m/s and 
the timespan is 4000 days. Radial velocity Fourier spectrum peak has amplitude K p and frequency n p . Fake associated planet has 
period T p = 2n/n p , semi-major axis a p , mass M p , and eccentricity e p . 









Theory 






Simulations 






Ex. 


Parameters 




r p (d) 


K p (m/s) 


T p (d) 


K p (m/s) 


a„ (AU) 


M p 


e„ 


1 


7=0° 


71 


8069 


- 


8032.8 ± 2.6 


6922.0 ± 0.3 


9.97 


1.05 M Q 


0.00 




e = e b = 


2n b — 3n 


222.6 


0.961 


223.2 ± 0.5 


0.96 ± 0.06 


0.720 


18.7 M e 


0.00 


2 


i = 30° 


77 


8069 


- 


8037.6 ± 2.7 


6916.1 ±0.3 


9.97 


1.05 M Q 


0.00 




e = e b = 


2 nj, — 3 n 


222.6 


0.837 


222.9 ± 0.6 


0.83 ± 0.06 


0.719 


16.0 M e 


0.00 


3 


i = 60° 


n 


8069 


- 


8048.3 ± 2.7 


6904.5 ± 0.3 


9.98 


1.05 M Q 


0.00 




e = e b = 


1n b — 3 n 


222.6 


0.541 


222.3 ± 0.9 


0.52 ± 0.06 


0.718 


10. 16 M e 


0.00 


4 


7=0° 


71 


8069 


- 


8035.1 + 2.8 


6958.1 ±0.4 


9.97 


1.05 Mq 


0.10 




e = 0.1 


2 M;, - 3 77 


222.6 


0.869 


222.6 ± 0.6 


0.90 ± 0.06 


0.719 


17.52 M ffl 


0.00 




e b = 0.1 


«i, — 3 n 


485.5 


0.585 


485.4 ± 3.7 


0.65 ± 0.07 


1.209 


16.01 M e 


0.00 


5 


7=0° 


n 


8069 


_ 


8026.8 ± 2.8 


6959.4 ± 0.4 


9.97 


1.05 M Q 


0.10 




e = 0.1 


n b - 3 n 


485.5 


1.149 


485.2 ± 1.8 


1.27 ±0.07 


1.208 


31.9 M e 


0.00 




e b = 0.2 


2 — 3 77 


222.6 


0.804 


222.2 ± 0.6 


0.84 ± 0.06 


0.718 


16.2 M e 


0.00 


6 


7=0° 


77 


8069 




8074 ± 6 


6993 ± 5 


10.02 


1.06 M Q 


0.10 




e = 0.1 


n b - 3n 


485.5 


1.673 


476.8 ± 1.4 


1.53 ±0.07 


1.194 


38.5 M e 


0.00 




e b = 0.3 


2 n b — 3 77 


222.6 


0.701 


222.6 ± 0.8 


0.68 ± 0.06 


0.719 


13.36 M ffi 


0.00 


7 


7=0° 


77 


8069 




8103 ±6 


7020 ± 6 


10.06 


1.07 Mq 


0.10 




e = 0.1 


7J(, - 3 77 


485.5 


2.135 


477.6 ± 1.1 


1.98 ±0.07 


1.196 


50.1 M e 


0.00 




e b = 0.4 


2 ftj — 3 77 


222.6 


0.569 


223.1 ±0.9 


0.57 ±0.06 


0.720 


11.2 M e 


0.00 


8 


i = 0° 


77 


8069 




7928 ± 2 


7550.1 ±0.4 


9.88 


1.05 M G 


0.39 




e = 0.4 


2 Mi — 6 71 


242.7 


1.104 


241.1 ±0.3 


2.60 ± 0.07 


0.758 


51.4 M e 


0.12 




e b = 0.05 


2 — 8 71 


258.2 


0.840 


253.7 ± 0.4 


2.52 ±0.07 


0.784 


51.1 M e 


0.00 



Table 2. Triple system of solar mass star and binary on circular coplanar orbits: parameters of planets mimicked by the binary. 



Ex. 


M, (M Q ) 


M 2 (M Q ) 


a (AU) 


a b (AU) 


T(y) 


K„ (m/s) 


t„ (y) 


a p (AU) 


M p (M e ) 


1 


0.3 


0.3 


10 


1.0 


10.26 


0.541 


0.685 


0.777 


16.00 


2 


1.0 


1.0 


10 


1.0 


10.26 


0.987 


0.375 


0.520 


23.90 


3 


1.0 


1.0 


50 


5.0 


114.7 


0.447 


4.197 


2.60 


23.90 


4 


1.0 


1.0 


50 


2.5 


114.7 


0.039 


1.482 


1.30 


1.494 


5 


1.0 


1.0 


100 


10.0 


324.4 


0.312 


11.86 


5.20 


23.90 


6 


1.0 


1.0 


100 


5.0 


324.4 


0.028 


4.197 


2.6 


1.494 



this system is currently only 60 m/s (Eggenber geret al1l2007l) . 
hence this effect is not detectable. 

4.4. Planets in binaries 

In|M orais & Correia! d2008l) . we reviewed all currently known 
planets of stars which are part of a moderately close binary sys- 
tem (a S5 100 AU). Our aim was to see whether these could be 
fake planets due to a binary composed of the companion star, 
Mi, and another unresolved star, M^. We computed the param- 
eters of the hidden binary c omponent (M? and ah) that could 
mimic these planets (Table l. lMorais & Correia d2008l) ) and saw 
that Mi would either be too massive to be realistic or at least too 
massive to be unr eported. 

The results in lMorais & Correial d2008l) are only valid in the 
context of triple star systems on circular coplanar orbits. Here, 
we saw that when the triple system's orbits are eccentric, the 
magnitude of the binary's effect can increase with respect to the 
circular case (Sees. 4.1 and 4.3). Therefore, in the eccentric case 
a planet can be mimicked by a less massive binary than in the 
circular case. 

The most interesting cases in Table 1 dM orais & Correial 
2008), namely yCep and HD 196885, are binary systems with 
eccentric orbits. The star yCep A orbit s yCep B with e = 0.4 1 
and / = 119° (Neuha user et al.1 12001 iMugrauer et all 120071) . 
while t abs x 772 dTorresll2007l) . The star HD 196885 A orbits 
HD 196885 B with e = 0.46 but / is unknown while t obs ~ T/4 



dCorreiaetalJ l2008). However, extrapolating from the results in 
Sect. 4.1, if there were hidden companion stars close to y Cep B 
or HD 196885 B, we would expect to be able to detect, in each 
case, several nearby peaks and not single peaks at the planets' 
orbital frequencies. 

Moreover, the stars yCep A+B have been directly imaged 
and y Cep B's brightness is consis tent with that of a single M4- 
type dwarf (Neuhau ser et al.l l2007). Additionally, we performed 
some 3-body fits to the data using a model with HD 196885 A+B 
and an unresolved star nearby HD 196885 B. We saw that these 
fits were worse (residuals at least 3 1 m/s and poorly constrained 
orbital elements) than fits using a model with the binary star sys- 
tem (HD 196885 A+B) and a planet com panion to HD 196885 A 
(residuals of 1 1 m/s in agreement with ICorreia et alJ (|2008)). 
Therefore, the planets around yCep A and HD 196885 A, re- 
spectively, are not likely to be due to hidden star companions to 
yCepB or HD 196885 B. 

5. Conclusion 

We have studied a triple system composed of a star and a binary 
system. The star and binary system have eccentric and inclined 
orbits. This is an extension of earlier work where we assumed 
that the star and binary system are on circular coplanar orbits 
dMorais & Correiall2008l) . 

We demonstrated that if we are unaware of the binary sys- 
tems presence (one or even both components may be unresolved 
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for instance because they are faint M stars) we may then be led 
to believe that the star has one or even two planet companions. 
Although the radial velocity variations due to a binary are dis- 
tinct from those due to planet(s), in practice, the measured effect 
depends on the instrument's precision and the observation time 
span. 

We have shown that because of the limited instrumental pre- 
cision, we may only detect periodic terms with well separated 
frequencies hence we mistake these for planet(s). We also saw 
that, if the observation time span is shorter than the star's long- 
period motion around the binary's center of mass, we may not 
be able to resolve terms with nearby frequencies, which means 
that we cannot distinguish fake planet(s) from a binary. 

We have also demonstrated that the binary's effect is more 
likely to be mistaken for planet(s) when the radial velocity os- 
cillations are composed of large dominant periodic terms with 
well separated frequencies. This is more likely to happen in 
the case of coplanar orbits observed equator-on. Moreover, we 
have seen that when the orbits are eccentric, the magnitude of 
the binary's effect can increase with respect to the circular case. 
Nevertheless, our model is valid for any triple system's configu- 
ration. 

We have presented an example of a binary system that af- 
fects a nearby star's motion. When the star's long-period motion 
has low to medium eccentricity, the binary can mimic planet(s) 
of 10 to 50 Earth masses, which are near the current detection 
limit. However, if the long-period motion has a high eccentric- 
ity we are more likely to detect multiple signals with very close 
frequencies and therefore reject the planet hypothesis. An excep- 
tion occurs when the long period motion has a high eccentricity 
but its period is poorly constrained (due to a short observational 
time span) in which case we may mistake the binary for a planet 
of about 100 Earth masses. 

We showed that when the binary has an eccentric orbit it 
can mimic two planets with periods approximately in the ratio 
2/1. Therefore we may be misled to think that we found plan- 
ets in the 2/1 orbital resonance when we have an eccentric hid- 
den binary. This is somehow analogous to the case studied by 
Anglada-Escude et al. (2010) where two planets on circular or- 
bits in the 2/1 orbital resonance can be mistaken for a single 
planet with an eccentric orbit. However, our scenario is probably 
more realistic since planets in orbital resonances are not likely 
to have circular orbits. 

We propose that new planet detections in close bi- 
nary systems, especially Earth-sized objects that are the tar- 
gets of the planned search program EXPRE SSO/CODEX 
dD'Odorico & COD EX/ESPRESS O Team! |200l . be checked 
carefully because they could indeed be artifacts caused by a hid- 
den binary. This could be done by comparing fits to the data us- 
ing (1) a model composed of the binary star system and a planet 
with (2) those obtained for a model composed of a hierarchical 
triple star system. If the fits of (2) are at least as good as (1), then 
the hidden binary hypothesis should be considered. 
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